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The MINOS long baseline experiment has been collecting neutrino beam data since March 
2005 and has accumulated 3 x 10 20 protons-on-target (POT) to date. MINOS uses Fermi- 
lab's NuMI neutrino beam which is measured by two steel-scintillator tracking calorimeters, 
one at Fermilab and the other 735 km downstream, in northern Minnesota. By observing 
the oscillatory structure in the neutrino energy spectrum, MINOS can precisely measure 
the neutrino oscillation parameters in the atmospheric sector. From analysis of the first 
year of data, corresponding to 1.27 x 10 20 POT, these parameters were determined to be 
|Am§ 2 | = 2.74±°;^ x 10" 3 eV 2 /c 4 and sin 2 (26> 2 3) > 0.87 (68% C.L.). MINOS is able to mea- 
sure the neutrino velocity by comparing the arrival times of the neutrino beam in its two 
detectors. Using a total of 473 Far Detector events, (v -c)/c = (5.1 ± 2.9) x 10 -5 (68% C.L.) 
was measured. In addition, we report recent progress in the analysis of neutral current events 
and give an outline of experimental goals for the future. 



1 Introduction 

It is now well established that neutrinos have non-zero masses and that neutrinos mix. Their 
weak interaction eigenstates (or "flavour" eigenstates) v a are related to their mass eigenstates 
Vi by a unitary transformation U: 



U is called the PMNS ^ matrix. Neutr inos are created and detected by weak interaction 
processes but their propagation in free space is described by their mass eigenstates, causing 
relative phases to change. This leads to the phenomenon of neutrino oscillations. 

MINOS is a long baseline neutrino oscillation search based at FNAL. The neutrino beam 
created in the NuMI beamline is sampled in two locations, ~lkm from the beam production 




(1) 



Correspondence: t .rauf erl@physics . ox. ac .uk 



Target Hall 

^ — 



Decay Pipe 



Muon Monitors 



Target 




Hadron Monitor 



Figure 1: The NuMI beamline components. Positively charged mesons produced in the target are focused by 
two magnetic horns and subsequently decay in a 675 m-long evacuated volume producing predominantly . The 
remaining hadrons are stopped in a beam absorber at the end of the decay pipe. 



target underground at Fermilab, and again 735 km downstream in the Soudan Underground 
Laboratory, in Minnesota. 

Using this setup, MINOS measures the oscillation parameters [Am^jand sin 2 (26>23)to world 
leading precision. In addition, MINOS searches for sub-dominant — > u e oscillations, oscilla- 
tions into sterile neutrinos, and v — > V transiti ons . The MINOS Far Detector can also be used 
to detect neutrinos created in the atmosphere. ^Sl 

2 Experimental setup 

2.1 The NuMI beamline 

A schematic of the NuMI neutrino beamline is shown in Figure [TJ NuMI uses protons with a 
momentum of 120 GeV extracted from Fermilab's Main Injector accelerator. They impinge a 
95.4 cm long segmented graphite target producing secondary particles, mainly ir and K mesons. 
The positively charged mesons are focused using a system of two pulsed, parabolic magnetic 
focusing elements, called "horns". The secondaries subsequently decay in a 675m long, 2m 
diameter evacuated decay volume producing neutrinos. 

Hadrons reaching the end of the decay volume without decaying are stopped in an beam 
absorber following the decay pipe. The beam absorber consists of a water-cooled aluminium core 
surrounded by a layer of steel blocks and an outer layer of concrete. The remaining muons from 
the meson decays are stopped in the ~300m of rock that separate the MINOS Near Detector 
hall from the beam absorber. 

The target position relative to the first horn and the horn current are variable. For most 
of the collection of data used in the oscillations analysis presented here, the target was inserted 
50.4 cm into the first horn to maximize neutrino production in the 1-3 GeV energy range. The 
data described here were recorded in this position, between May 2005 and February 2006, and 
correspond to a total of 1.27xl0 20 POT. The charged current (CC) neutrino event yields at the 
ND are predicted to be 92.9% 5.8% F M , 1.2% u e and 0.1% V e . 

2.2 The MINOS detectors 

The MINOS detectors are designed to be as similar as possible while operating in very different 
conditions. They are steel-scintillator sampling calorimeters, magnetized to ~1.3T allowing to 
measure particle momentum using track curvature as well as range. The scintillator planes are 
divided into 4. 1 cm wide strips rotated by 90° on subsequent planes to enable 3-dimensional 
event reconstruction. The scintillator light is captured by embedded wavelength shifting (WLS) 



fibres and transported to the edge of the detectors, where the optical signal is converted using 
Hamamatsu (Near Detector) and M16© (Far Detector) photomultiplier tubes. 

Below 10 GeV, the hadronic energy resolution was measured to be 56%/ \JE [GeV] ©2% and 
the EM resolution was measured to be 21A%/ y/ E[GeV] © 4.1%/ E[GeV] ™ The muon energy 
resolution AE^/En varies smoothly from 6% for above 1 GeV where most tracks are contained 
and measured by range, to 13% at high energies, where the curvature measurement is primarily 
used. 



3 Monte Carlo tuning 

MINOS took data with 6 beam configurations obtained by varying the target position and 
the current in the focusing horns. The Monte Carlo simulations of neutrino fluxes strongly 
depend on the underlying models of hadron production in the target, which are presently poorly 
constrained at MINOS energies. The nominal simulations, based on FLUKA05, yielded energy 
spectra which did not match the high-statistics data in the Near Detector. A better agreement 
was achieved by smoothly adjusting the p z and pr of hadrons (mostly pions) produced off the 
graphite target. The resulting spectra, which describe the data much more closely than the 
nominal simulations, are shown in Figure El 
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Figure 2: Data and Monte Carlo energy spectra and their ratios are shown for six deferent beam configurations. 
The blue line corresponds to the untuned, the red line to the tuned hadron production model. The tuning 

significantly improves data/MC agreement in all cases. 



4 Uu disappearance analysis 



4-1 Event selection and extrapolation 

In the Far Detector (FD), the signal from eight scintillator strips is read out by the same PMT 
pixel. Therefore, the initial step in the reconstruction of the FD data is the removal of the 
eightfold hit-to-strip ambiguity using information from both strip ends. In the ND, timing and 
spatial information is first used to separate individual neutrino interactions from the same spill. 
Subsequent reconstruction is done in the same way in both detectors. Tracks are found and 
fitted, and showers are reconstructed to be combined to events. For v„ CC events, the total 
reconstructed event energy is obtained by summing the muon energy and the visible energy of 
the hadronic system. To prevent human biases when assessing the oscillation analysis results, a 
blinding mechanism was applied to the FD data set. This procedure hid a substantial fraction 
of the FD events with the precise fraction and energy spectrum of the hidden sample unknown. 
Events are pre-selected in both detectors, by requiring total reconstructed energy below 30 GeV 
and a negatively charged track. The track vertex must be within a fiducial volume such that 
cosmic rays are rejected and the hadronic energy of the event is contained within the volume of 
the detector. The pre-selected Un event sample is predominantly CC with a 8.6% neutral current 
(NC) event background estimated from Monte Carlo (MC) simulations. The fiducial mass of 
the FD and ND is 72.9% and 4.5% of the total detector mass respectively. 

A particle identification parameter (PID) incorporating probability density functions for the 
event length, the fraction of energy contained in the track and the average track pulse height 
per plane provides separation of Un CC and NC events. The PID is shown in Figure[3]for ND 
and FD data overlaid with simulations of NC and CC events. Events with PID above -0.2 (FD) 
and -0.1 (ND) are selected as being predominantly CC in origin. These values were optimized 
for both detectors such that the resulting purity of each sample is about 98%. The efficiencies 
for selecting CC events in the fiducial volume with energy below 30 GeV are 74% (FD) and 
67% (ND). 
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Figure 3: Data and tuned MC predictions for the PID variable in the ND (top) and FD (bottom). The arrows 
depict the positions of the selection cuts. The FD MC distribution for CC events uses the best fit parameters 

discussed in the text. 

The measurement of the energy spectrum at the ND is used to predict the unoscillated 
spectrum at the FD. Fits to the ND data yield tuning parameters for the predicted neutrino 
flux. These fits are based on parameterisations of the secondary pion production at the NuMI 
target as a function of xf and px as described in Sectional The FD prediction must also take into 
account the ND and FD spectral differences that are present, even in the absence of oscillations, 




due to pion decay kinematics and beamline geometry. This is achieved using the Beam Matrix 
method. It utilizes the beam simulation to derive a transfer matrix that relates the neutrinos 
in the two detectors via their parent hadrons. The ND reconstructed event energy spectrum is 
translated into a flux by first correcting for the simulated ND acceptance and then dividing by 
the calculated cross-sections for each energy bin. This flux is multiplied by the transfer matrix 
to yield the predicted, unoscillated FD flux. After an inverse correction for cross-section and FD 
acceptance, the predicted FD visible energy spectrum is obtained. The oscillation hypotheses 
are then tested relative to this prediction. A distinct extrapolation method, referred to as ND 
Fit was also applied to the data, yielding similar results. 

In total, 215 events are observed below 30GeV compared to 336.0±18.3(stat.)±14.4(syst.) 
events expected in the absence of oscillations. The systematic error is most relevantly due to 
NC contamination, ND to FD normalization and the hadronic shower energy scale. In the 
region below lOGeV, 122 events are observed compared to the expectation of 238.7±15.4±10.7. 
The observed energy spectrum is shown along with the predicted spectra for both extrapolation 
methods in Figure HI 
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Figure 4: Comparison of the Far Detector spectrum with predictions for no oscillations for both analysis methods 
and for oscillations with the best fit parameters from the Beam Matrix extrapolation method (left canvas). The 
estimated NC background is also shown. The last energy bin contains events between 18-30 GeV. The right canvas 
shows the ratio of data and best fit over the unoscillated predictions. 



4-2 Oscillation analysis 

Under the assumption that the observed deficit is due to — * v T oscillations, a x 2 fit is 
performed to the parameters |Am 2 2 | and sin 2 (2^23) using the expression for the survival 
probability: 

P(u, -„„) = !- sin 2 (2c? 23 ) sin 2 (^f^) (2) 

where L is the distance from the target, E is the neutrino energy, and jAm^l is the atmospheric 
mass splitting. The fit included the systematic uncertainties mentioned above as nuisance pa- 
rameters as well as the small contribution from selected v T events produced in the oscillation 
process. The resulting 68% and 90% confidence intervals are shown in Figure [5] as determined 
from A% 2 =2.3 and 4.6, respectively. The best fit parameter values are: 

|Am 2 2 | = (2.74i° -i 4 6 )xl0- 3 eV 2 /c 4 (3) 

and 

sin 2 (2#2 3 ) > 0.87 (4) 

at 68% C.L. with a fit probability of 8.9%. At 90% C.L.(2.3K|Am§ 2 |< 3.43) x 10" 3 eV 2 /c 4 , 
and sin 2 (26>23) > 0.78. The data and best fit MC are shown in Figure HI 
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Figure 5: Confidence intervals for the fit using the Beam Matrix method including systematic errors. Also shown 
are the contours from the previous highest precision experiments. IS12UQI 



If the fit is not constrained to be within the physical region, the best fit is at |Am3 2 |=2.72 x 
l(T 3 eV 2 /c 4 and sin 2 (26>2 3 ) = 1.01, with a decrease in X 2 of 0.2. 

It is expected that the systematic uncertainties will be reduced with additional data. More 
details of this analysis are available inLUJ An update on this result using 2.58 xlO 20 POT is 
expected during the Summer 2007. 



5 Neutrino time-of-flight analysis 

MINOS uses GPS synchronized clocks to timestamp neutrino interactions in both detectors. 
This enables the measurement of the neutrino time-of-flight over a distance of 734 km and thus 
the determination of the neutrino velocity. Similar terrestrial experiments performed in the past 
used much shorter baselines of ~ 500 m and higher beam energies (> 30 GeV). 

The time of each PMT hit is recorded by the detector's clock to a precision of 18.8 ns in 
the Near Detector and 1.6 ns in the Far Detector. The time of the earliest hit of each event 
is taken to be the time of the neutrino interaction. The interaction times t\ and ti in the two 
detectors are corrected for known offsets and delays, which were determined using test stand 
measurements. In addition, the time of the beam extraction signal is subtracted from both 
times. The beam extraction signal has a fixed relation to the arrival of neutrinos in the MINOS 
detectors. All times are therefore measured relative to this reference. 

The NuMI beam pulse is not instantaneous, but has a duration of 9.7 \xs with a five-batch 
or six-batch intensity profile depending on the accelerator running mode. The Near Detectors 
measures this time-intensity profile with neutrino interactions to high precision. The measured 
time-intensity profile forms a probability density function which is folded with a Gaussian dis- 
tribution with a width of a = 150 ns to account for the uncorrelated jitter of the two GPS 
clocks. The resulting distribution P{t) describes the predicted arrival time distribution at the 
Far Detector (shown as a solid line in Figure [6] for the five- and six-batch modes separately) . 

For the time-of-flight measurement, 473 neutrino- induced events in the Far Detector were 
used. The time of each event was compared to the predicted arrival time distribution. The 
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Figure 6: Time distribution of FD events relative to prediction after fitting the time-of-fiight. The top plot shows 
events in 5-batch spills, the bottom 6-batch spills. The solid lines show the normalized prediction curves. 

time-of-flight r was found by maximizing an unbinned log-likelihood function: 

L = 5>P(4~t). (5) 

i 

The distribution of measured event times together with the predicted distribution for the best 
fit r is shown in Figure [6j The time-of-flight of neutrinos was measured to be 

2449.223 ± 0.032(stat.) ± 0.064(syst.) /js 68% C.L. (6) 

Comparing to the MINOS baseline this translates to 

= 5-1 ± 2.9 (stat. + sys.) x 1(T 5 68% C.L. (7) 
The systematic error is due to uncertainties on the timing delays and offsets. 

6 Progress in Neutral Current analyses 

Analyses of neutral current neutrino interactions in MINOS are currently in progress. Neutral 
current interactions are interesting for several reasons. They form an important background 
to the — > v T oscillation analysis described in these proceedings and their cross-sections at 
the energies relevant to MINOS are not very well known. Furthermore, an observed deficit of 
neutral current events at the Far Detector could be evidence for light sterile neutrinos. 

Neutral current interactions are selected using three event quantities: the event length, the 
number of tracks in the event and the track extension. In events where both tracks and showers 
were found, the track extension measures how much longer a reconstructed track is compared 
to the hadronic shower it is accompanied by. Distributions of the selection variables in Near 
Detector data and Monte Carlo are shown in Figure [7] (left canvas). 

Using these variables, an energy spectrum of neutral-current-like events was produced. This 
is shown in the right canvas of Figure[7l The Monte Carlo is shown as a red line with a systematic 
error band; the background due to wrongly selected charged current events is shown as a blue 
hatched distribution. Within the estimated systematic errors, due to flux, cross-section and 
energy scale uncertainties, data and Monte Carlo agree well. 

The extrapolation of these Near Detector results to the Far Detector and a fit for oscillations 
to sterile neutrinos are currently being worked on. Results from this analysis are expected later 
this year. 
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Figure 7: Neutral current event selection variables in the Near Detector (left). The data is shown as black markers, 
the Monte Carlo in red. The blue distribution shows the charged current background. The right canvas shows 
the reconstructed energy spectrum for selected neutral current events. The Monte Carlo (red) is shown with 

systematic errors. 



7 Future Prospects 

In addition to the analyses reported here, MINOS is pursuing a v„ — > v e oscillation analysis 
in order to measure or constrain the as yet unknown mixing angle #13. With its expected 
final statistics, MINOS will potentially be able to improve on the current best limit from the 
CHOOZ^ 

experiment. 

Other areas of interest include appearance and disappearance measurements of anti-neutrinos 
in the Far Detector as well as several non-oscillation analyses using the large number of neutrino 
interactions in the Near Detector. 
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